An Anderson Transition of the Plasma Oscillations of ID Disordered Wigner Lattices 



OO 

o 
o 

(N 

C 



I 

CZ3 



cd 



T3 

C 

o 

> 

o 

o 
\o 
o 

-I— > 

c3 



c 
o 
o 



X 



Shimul Akhanjeqj and Joseph Rudnidtj] 
Department of Physics, UCLA, Box 951547 Los Angeles, CA 90095-1547 

(Dated: June 2, 2008) 

We report the existence of a localization-delocalization transition in the classical plasma modes 
of a one dimensional Wigner Crystal with a white noise potential environment at T — 0. Finite 
size scaling analysis reveals a divergence of the localization length at a critical eigenfrequency. 
Further scaling analysis indicates power law behavior of the critical frequency in terms of the relative 
interaction strength of the charges. A heuristic argument for this scaling behavior is consistent with 
the numerical results. Additionally, we explore a particular realization of random-bond disorder in 
a one dimensional Wigner lattice in which all of the collective modes are observed to be localized. 
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Recently, there has been considerable interest in 
the localization phenomena of disordered, one dimen- 
sional(lD), electronic tight-binding models. Moreover, 
it has been established that an Anderson transition can 
occur in lower dimensional systems containing nearest- 
neighbor interactions and random on site energies ex- 
hibiting long range correlations in its power spectrum [l| . 
Relatively little is firmly understood about the classical 
oscillator models with long range interactions and disor- 
der. Wigner Crystal(WC) [2] phases are one such exam- 
ple, having collective excitations that posess a different 
character from the weakly interacting regime, especially 
in the presence of disorder. Various authors have stud- 
ied plasmon localization in this context. However those 
models do not consider the true long range nature of 
the coloumb potential as discussed in this article [3j. In 
this extreme limit of electron-electron interactions, the 
dominant behavior of the collective modes can be treated 
classically due to the negligible overlap of the electronic 
wavefunctions at lower densities [4|]. 

In this paper we report on the existence of a 
localization-delocalization transition in the plasma os- 
cillations of a disordered ID WC at zero temperature. 
We provide strong evidence for the existence of such a 
transition based on numerical studies of both the local- 
ization length £ and the inverse participation ratio P. 
In particular we show that the transition occurs in a ID 
system with uncorrelated disorder in contrast to the mo- 
bility edge observed by Izrailev et al [l(. Data collapse 
at various system sizes confirms both scaling at the tran- 
sition and a critical exponent v associated with the diver- 
gence of the localization length. Additionally, we recover 
a scaling exponent of the critical eigenfrequency in terms 
of the relative interaction strength of the charges. We 
also explore a disordered model with randomness in the 
effective bonds between the charges. This can be viewed 
as an extension of Dyson's disordered chain problem [5| . 
In the case of this system no Anderson-like transition 
exists. Rather, all plasma oscillation eigenstates are lo- 
calized. In this regard, the system mirrors the behavior 
of the one-dimensional system with short-range disorder. 



Unlike that system, the localization length is not well de- 
scribed by the Thouless formula which follows from the 
mathematical simplicity of a dynamical matrix having 
tri-diagonal structure [6(. 

The first system under investigation, which we term 
Model A, consists of like charges placed in an external 
random potential with the charge neutrality condition en- 
forced by smearing a positive background over the chain. 
The Hamiltonian for Model A is given by: 
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We assume the low density limit in which Coloumb in- 
teractions dominate the kinetic energy term, resulting in 
a Wigner solidification of the particle array. V(x) is a 
white noise potential with the following Fourier decom- 
position, with N = L/4 Fourier components. 

N 

V(x) — A 2_] an cos(2nnx) + b n sin(27rnx) (2) 



The parameters J and A in Eqs. |1| and © are dimen- 
sionless coupling constants, which we utilize to define the 
dimensionless interaction strength k = J /A. The random 
variables a n and b n are chosen from different Gaussian 
distributions, having an independence of n, yielding a 
white noise power spectrum with a mean of /i = and 
a variance of a = 1. In Model B the external potential, 
V(x) = and the charges {Qi} are random variables 
with a distribution function P[Q], given by: 



P[Q] = <j Wa-m 



Wx < Q < W 2 
otherwise 



(3) 



The formalism used for computing the eigenfunctions 
of the collective modes involves the construction of the 
dynamical matrix D(R) . For a finite sized chain of length 
L, D(R) is an I x I symmetric matrix. In the harmonic 
approximation the full electrostatic potential is expanded 
in a Taylor series with the terms proportional to the sec- 
ond derivatives retained. The dynamical matrix has the 



structure 
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(4) 
where <p(x) is defined as the electrostatic potential be- 
tween two charges in a periodic image and the R's are the 
equilibrium positions of the charges to be determined by 
a numerical minimization. The dynamical matrix enters 
into the the eigenvalue equation for plasma eigenmodes 
as follows [7(: 



m e LU 2 u(R) + ^ D(R - R')u(R') = 
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where the functions u(R) are the lattice displacements 
from equilibrium. Periodic boundary conditions were em- 
ployed to eliminate surface effects, and lengths are scaled 
so that the size of the region occupied by the charges is 
unity. We make use of a closed form expression for the 
electrostatic potential, 4>{x), that possesses the correct 
periodicity, short distance behavior and that properly 
reproduces the standard Ewald potential result to high 
degree of accuracy. The explicit form that we utilized is 



4>{x) — 7r| csc(7ra;)| 



(6) 



After assigning charge values to the particles in an 
ordered array, the system must relax so that the total 
force on each particle is zero. We employed the Newton- 
Raphson method to numerically determine the particle 
coordinates at equilibrium. 

Given the relaxed array of charges, we are in a po- 
sition to calculate the dynamical matrix D(R) via (jil) 
and the solutions of the eigenfunction equation ([5| . The 
eigenfunction width w is a direct measure of the local- 
ization length £, taken as the mean squared deviation of 
the mimimum distance between lattice sites. A natural 
choice is: 
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(7) 



where Ui(k) is the amplitude of the ith eigenfunction at 
the site k. The one difficulty with this definition of the 
width is that a function localized to the boundaries of the 
interval would be described as extended. To correctly 
take into account the periodicity of the system and to 
avoid the above pitfall, we replace (k — I) 2 in the first 
line of (O by sin 2 7r(fc - 1)/tt 2 . 

Another important measure of spatial localization is 
the inverse participation ratio (IPR) P [8] , defined as 




FIG. 1: Typical extended and localized eigenfunctions(Model 
A, L=512) 



p = lJ2 mj)i 4 /(E i«iC?)i 2 ) a = lJ2 kooi 4 ( 8 ) 

The quantity P is a useful measure of which fraction of 
the oscillators are localized. As the thermodynamic limit 
is taken, P approaches zero for extended modes, while for 
localized modes the IPR can be as high as P — L. For 
both models discussed in this article, one must account 
for the numerical and statistical uncertainties on finite 
length chains that scale like oc l/y/L. Therefore, w and 

P were ensemble averaged over 50 realizations of disorder 
to smooth out these fluctuations. 

Typical eigenfunctions in the extended and localized 
regimes of Model A are shown in Fig. [1] Due to the 
limitation of performing diagonalizations of finite sized 
matrices, one can interpret an extended eigenfunction as 
having a width value, Wi = £ = L. The existence of 
criticality is indicated by a common crossing point for 
the values of w and P that are normalized at various 
system sizes. A genuine Anderson transition is usually 
accompanied by a diverging localization length at a par- 
ticular eigenvalue [9( and a similiar abrupt shift in the 
IPR at this same eigenvalue. For electronic tight binding 
systems, the critical eigenvalue that separates conducting 
and insulating phases is known as the mobility edge [l0( . 
Likewise, in Model A the plasma oscillations display two 
frequency regimes, one corresponding to localized and the 
other to extended lattice waves analogous to the mobil- 
ity edge. However the DC conductivity for charge trans- 
port is equal to zero as Wigner crystals are known to be 
Mott Insulators |2|. Furthermore, the interaction ma- 
trix elements of classical oscillator systems are not con- 
strained to even integer sites as usually implemented in 
tight binding systems with random hopping in order to 
preserve the particle-hole symmetry of the Hamiltonian. 
Therefore along with our required structural relaxation 
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FIG. 2: Inverse participation ratio vs. eigenvalues normalized 
at different system sizes showing transition point in Model A. 



procedure, the Models A and B fundamentally differ from 
those investigated by Bhatt and Zhou [ill ]. 

The IPR data is presented in Fig. [2] For all cases 
in Model A, the interaction strength parameter, k is set 
equal to 0.1. Note that for a range of low frequencies, the 
localization lengths saturate at £ = L, while at higher 
frequencies, £ — » 0. In fact, Fig. [2] Model A appears to 
be consistent with a transition between frequency regimes 
in which the lattice excitations shift between the general 
forms illustrated in Fig. [T] Additionally, Fig. [2]exhibits a 
common crossing point in Model A that is consistent with 
asymptotically diverging IPR values above the crossing 
point as a function of L, signifying localization. 

In another test for a transition, we attempt the appli- 
cation of finite size scaling analysis to the dependence of 
the localization length of the lattice modes of the Model 
A system 12|, |l3[ . The results of this analysis are shown 
in Fig. [3l We test for a scaling form for the localization 
length having the following dependence on key quantities 
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The function F a is a universal function present in phase 
transitions belonging to the same universality class, the 
precise nature of which has not yet been determined. One 
sees a clear collapse of the data at 1/v ss 0.45 ± 0.05. 
We performed a statistical analysis of the curves at dif- 
ferent system sizes to provide a quantitative measure of 




FIG. 3: Normalized localization length vs. scaled eigenvalues 
for different system sizes, showing data collapse at \jv ~ 
0.45 ± 0.05 (Model A) 



the "goodness" of the collapse. Our results yield a chi 
squared p value of p = 0.37, indicating no statistically 
significant differences between the collapsed curves. For 
fixed k, the critical frequency lo c is asymptotically in- 
dependent of system size in the limit of a large system. 
The dependence of the square of the critical frequency 
uj\ on the relative interaction strength was examined be- 
tween the values of n — 0.075 and k = 0.1, for the values 
of ui 2 c within the center of the band. In this regime the 
Coulomb force is close to the random force in strength. 
We performed a scaling analysis in this non-perturbative 
regime, assuming a scaling form of u) 2 — An s . We have 
determined 8 « 0.55 ± 0.008 from a linear fit with an R 
squared value of 0.999. 

We can attempt to qualitatively understand this ex- 
ponent by considering the Hamiltonian containing the 
random potential and the unscreened coloumb interac- 
tions. We speculate that the localization-delocalization 
can be attributed to the threshold at which the range of 
the Coloumb force is comparable to the length scale of 
the net disorder force. We define this length scale To. 
The explicit scale dependence of the Coloumb force and 
the pinning forces are given as 
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(10) 



(11) 



Clearly, the net random force should follow a random 
walk that scales with the number of pinning centers or 
in terms of the characteristic length scale, the length be- 
tween the particles. It follows that at the threshold, 
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(12) 



The critical eigenvalues will then scale as, 



- ~ K 3 / 5 



(13) 



We find it rather interesting that this value of S = 0.6 
compares well to our numerically determined value of 
S sa 0.55 ± 0.008. It is possible that this scale depen- 
dent behavior can be connected to the collective pinning 
of elastic manifolds at the Larkin scale 14 1. The low 
frequency dynamics of pinned elastic systems has been 
studied by Fogler, where it is argued that the localiza- 
tion length for plasmons is related to the Larkin length 
defined as the size R c of domains individually pinned by 
spatially random potential wells [3j. 

On the other hand, model B appears to admit of no 
Anderson transition as shown in Figs. [2] All oscillator 
systems contain the extended eigenmode corresponding 
to the w 2 = eigenvalue. For Model B as the thermody- 
namic limit is approached, £ asymptotically approaches 
this one extended mode while the rest of the modes are 
localized. Additionally, no crossing point exists in the 
IPR data, confirming an absence of criticality. We have 
checked this result for various distribution widths and 
mean values defined in ([3]), and no criticality has been 
observed. For both models one can attempt to under- 
stand the precise behavior of the eigenwidths in terms 
of the Lyapunov exponent 7 = l/2£ derived from the 
residue of the single particle Green's function [6(] . 
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In Eq. (EH, g{^ 2 ) is the density of normal modes 
and Cof is the matrix cofactor of the inverse resolvent. 
For systems with nearest-neighbor interactions the sec- 
ond term in Eq. (|14[) does not depend on the eigen- 
values and averages to a constant yielding the Thouless 
formula [6j. However, for systems with long-ranged in- 
teractions this term exhibits a nontrivial eigenvalue de- 
pendence. An Anderson transition would correspond to 
7(w 2 ) = 0. Ultimately, the existence of unscreened \/R 
interactions is simply not the sole factor responsible for 
the transition. Rather it is the complex and subtle inter- 
play of disorder and interactions. 

We conclude by noting that the localization properties 
of lattice waves has important consequences with regard 
to the energy and heat transport of real physical sys- 
tems. Recent articles have discussed a new class of ma- 
terials that contain structurally ordered charge configu- 
rations such as the copper-oxide chains of compounds like 
Nai+ x Cu02 [15j. These materials contain magnetic and 
thermodynamic properties that are consistent with some 
type of ID Wigner lattice formation. A possible thermo- 
dynamic signature of an Anderson transition observed in 



Model A would be present in the heat transport, given 
that a major source of the carriers of energy would arise 
from the collective excitations of the electronic lattice. 


Finally, we note that eigenfunction localization has 
been observed in random banded matrices with power 
law modulation by Mirlin et. al. [l7| and by Levitov 
[18j | in one-dimensional tight-binding Hamiltonians with 
long-range hopping. The Hamiltonian governing the dy- 
namics of plasmons in our model shares important fea- 
tures with those systems in that in all three cases the 
issue is the spatial structure of the eigenfunctions of ma- 
trices that incorporate randomness and that contain en- 
tries that decay as a power law in the difference between 
indices. However, key power laws are different, in that 
we observe an effective mobility edge in eigenvectors of 
an operator in which off-diagonal terms decay effectively 
as \i-j\ 



the case of 17|, |1 



rather than as the inverse first power, as in 
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